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Abstract 

Tunneling measurements of the temperature dependence of the electronic density of states (DOS) 
of ultrathin bilayers of Pb and Ag reveal that their superconducting energy gap, A(T), evolves 
similarly to BCS predictions despite the presence of a large anomalous DOS at the Fermi energy 
that persists as T — > 0. The gap ratio, , systematically decreases more than 20% below the 
BCS universal value of 3.52 in bilayers with the lowest superconducting transition temperatures, 
T c ; a behavior we deem super weakly coupled. A semi-quantitative model suggests that the reduced 
gap ratios result from a systematic depletion of the DOS available for pairing that occurs with the 
growth of the anomalous DOS at the Fermi energy. 
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The BCS theory of superconductivity predicts that the gap ratio, r\ = 2A /&:bT c , assumes 
a universal value of 3.52 Here, T c and A are the superconducting transition temperature 
and zero temperature energy gap, respectively. Experiments on a large number of elemen- 
tal superconductors show 77 to range from a little below 3.5 to as high as 5 j^, in rough 
accord with this prediction. Extensions of the BCS theory into the strong electron phonon 
coupling regime |3j reveal that rj reflects the strength of the microscopic coupling constant 

n 

and show that stronger coupling correlates with larger 77 |2j. This correlation is often used 
to classify novel superconductors as weakly or strongly coupled even when the microscopic 
electron pairing interactions responsible for their superconductivity are unknown jj]. Here, 
we present a superconducting system, ultrathin bilayers of superconductor (S) and normal 
metal (N), in which 7] can be driven well below 3.5. We deem this surprising behavior "super 
weakly" coupled superconductivity. 

Ultrathin SN bilayers, are expected to behave similarly to pure superconducting layers 
according to quasi- L** p r0 ; imity effect tKeonesB . Their gap ratio should conform to 
i] > 3.5 and their Density Of States (DOS) should exhibit an energy gap devoid of states, 
centered on the Fermi energy, Ep, provided they are thinner than the superconducting coher- 
ence length, £o- The presence of the N layer simply decreases the bilayer's T c |6j. Tunneling 
experiments, however, indicate that bilayers with ultrathin S layers are not conventional. 
They possess an anomalous DOS within the energy gap region which becomes more pro- 
nounced with increasing normal layer thickness, 0,0,0]. In order to fit the whole DOS 
it is necessary to presume the coexistence of two distinct populations of quasiparticles, one 
giving rise to a BCS-like DOS and the second contributing the subgap states ~\. 

Here we show that the temperature dependence of the average energy gap, A(T), closely 
resembles the BCS form indicating that the growth of the superconducting condensate is 
conventional and uniform through the bilayers. In the bilayers with the lowest T c 's, however, 
A(T) asymptotes to an unexpectedly low value so that rj < 3.52. Concomitantly, the 
anomalous DOS at the Fermi energy appears to persist in the T —>■ limit and is largest 
in bilayers with the lowest rj. Present theories indicate that these states correspond to 



HAH. By assuming 



quasiparticles that are uncoupled from the superconducting phase [101. Ill 
that their decoupling reduces the DOS available for pairing it is possible to account for the 
anomalously low values of 77 using simple formulae for T c and A and thus, the super weakly 
coupled behavior. 
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The Pb/Ag bilayer films employed here were quench condensed from vapor onto T = 8 K 
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glass substrates in the UHV environment of a dilution refrigerator based cryostat 
Prior to cryostat mounting a substrate was cleaned, fire polished and patterned with Au/Ge 
contacts and Al (Alloy 2024) strips that served as the tunnel junction counterelectrodes. The 
quench condensed bilayers were formed by depositing an electrically discontinuous Pb film 
first and then a series of Ag layers. The mass per unit area of each deposition was measured 
with a quartz crystal microbalance and converted to a thickness using the bulk density of 
each element. After each Ag deposition, the tunnel junction conductance was measured as 
a function of temperature using standard 4-terminal, low frequency AC techniques. In this 
way, a series of Ag/Pb bilayers with fixed Pb thickness dpi, and increasing Ag thickness 
d-Ag were fabricated and measured in situ with the same tunnel junction barrier, without 
breaking vacuum or warming. This paper presents data on a series of 5 bilayers with dpi, 
= 4.0 nm, 6.7 nm < d&g < 19.3 nm, and 0.72 K < T c < 2.55 K. This series exhibits 
characteristics consistent with other series' of bilayers with 1.4 nm < dpi, < 4.0 nm. 

The normalized tunneling DOS of a bilayer, N$(E), is related to the normalized conduc- 
tance of the tunnel junctions through [lfij ]: 

df(E + eV) 
d(eV) 

where E is the energy measured relative to the Fermi energy, Ep, N$(E) = N$(E) / Nn(0) is 
the superconducting DOS normalized by the normal state DOS at Ep. df(E + eV)/d(eV) 
is the derivative of Fermi function and V is the voltage across the junction. 

The temperature dependence of Gj, shown for a T c = 2.0 K bilayer in Figs, la and 
lb, primarily follows BCS predictions with one very significant deviation. In accord with 
BCS, symmetric peaks grow and a depression between the peaks deepens as T decreases 
below T c . Most importantly, a well defined gap edge at « ±0.29 meV develops at the lowest 
emperatures. Contrary to BCS, the depression does not drop to zero within the gapedges 
1^ | . Rather, it develops into a symmetric, linear cusp with a finite intercept (see Fig. lb). 
Experiments on other bilayers and trilayers assure us that this subgap conductance reflects 
an intrinsic feature of the bilayer DOS and not junction leakage Q|. Moreover, investigations 



°^ T) = L Ns{E) d(eV) dE 



by other groups on other mesoscop, hybrid structures such as thin Au/Nb bilayers^ and 
inhomogeneous TiN films reveal a similarly anomalous subgap DOS. 

The behavior above suggests that the finite DOS at Ep persists to arbitrarily low temper- 
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FIG. 1: (a) Gj vs. V of a bilayer film (d Pb /d Ag = 4.0 nm/9.1 nm) at T = 0.06, 0.1, 0.3, 0.6, 1.0, 
1.4, 1.9 K. (b) Data and fits in the subgap region, (c) The inferred normalized energy gap A/Ao 
vs. T/T c . 

atures. To check this further, we measured the zero voltage bias conductance, Gj(0, T), as a 
function of temperature. Gj(0, T) is proportional to the DOS within ~4 kT of E F . Data for 
five of the bilayers with 0.72 K < T c < 2.55 K, are compared to the BCS dependence in Fig. 
2a. With decreasing T, Gj(0, T) shows an abrupt slope change at T c signalling the opening 
of the energy gap and subsequently, decreases monotonically in agreement with BCS. At 
the lowest temperatures, however, Gj(0,T) is higher than the BCS prediction and decreases 
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FIG. 2: (a) Semi-logarithmic plot of Gj(0) vs. T/T c for Pb/Ag bilayers with dpb = 4 nm and 
d-Ag = 6.7, 9.1, 12.4, 15.6, 19.3 nm. The dashed line is the mean field theoretical prediction of 
Gj(0) vs. T. (b) Same data as in (a) with fits (dashed lines). (c)Energy gap A vs. T obtained 
from the fits in b) (dpb = 4.0 nm and d^g = 6.7 nm (solid triangle), 9.1 nm (solid square), 12.4 
nm (empty circle), 15.6 nm (empty square) and 19.3 nm (empty triangle)). 

nearly linearly with temperature rather than exponentially. In fact, the Gj(0,T) appear to 
asymptotically approach finite values. The limiting values, Gj(0, 0), are larger for the lower 
T c bilayers. 

Both Gj(V, T) and Gj(0, T) can be fit by introducing temperature dependence into a two 
component form for the DOS that was introduced previously to fit tunneling data at a single 
low temperature T =60 mK «T C jjj. This DOS consists of a broadened BCS DOS at high 
energies and a linear DOS with an offset at energies within the gap. The former corresponds 
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to BCS-like quasiparticles and t 
localized in the normal regions 



re latter corresponds to quasiparticle states that are quasi- 
Ill [ill Q| • Those contributing to the linear portion (aE) 



are partially localized while those giving the offset (/3) are completely localized within the 
N layer |v|j . Explicitly, 

_ f aE + (3 ■ E < E c 

N S (E) = { _ 

where E c is the crossing energy separating the two DOS forms. N$(E) is a superposition of 
BCS DOS' Ng CS (E, A) = Re(E/\ / E 2 — A 2 ) using a log-normal distribution of A Q: 

1 r*SrB08,n . , 



NI(E)=^=—-[ Ni cs (E,A)exp( 
\/2ir koVo 



^(ka)h b v ' ' ^ v 2(£;a) 2 y A 

. The product fccr characterizes the width of the log-normal gap distribution and A is the 
most probable A. The fits to the data shown in Fig. lb were optimized to capture the 
gap edges. They matched the peaks to within 10%. The inability to capture both features 
simultaneously probably results from inaccuracies in the log normal approximation for the 
gap distribution. A was adjusted for each temperature, a, (3 and (ka) were presumed 
temperature independent and held to the values that produced the best fit at T =0.06 K. 
We are confident that this presumption is reasonable because of the physical interpretation 
of the parameters Q|. It gives high quality fits to the lowest temperature data where these 
parameters exert their greatest influence implying that a subgap DOS persists in the T — ► 
limit. As shown in Fig. 1(c), A(T) agrees well with the BCS form (the solid line) and gives 
an rj = 3.57, which is close to the BCS value (3.52). 

Fits to the zero bias conductances for all of the bilayers, yield A(T) and a good 
estimate of the density of electronic states at Ep that exists in the zero tempera- 
ture limit. The fitted curves shown in Fig. 2b were calculated using: Gj(0) = 
fl^ N$(E)[df(E + eV) / d(eV)]v=odE , presuming that a, (3 and (ka) are constant (see ear- 
lier discussion) and adjusting A at each temperature point along a curve. At low temper- 
atures, A must approach a constant value, Ao, in order for the fits to reproduce the low 
temperature linear dependence of Gj(0). These Ao agree with the most probable energy gap 
A obtained from fits to Gj(V) at 60 mK to within ~15% in the worst case. The resulting 
A(T) for all of the bilayers, which are shown in Fig. 2c resemble the BCS form. Thus, the 
fits to Gj(V,T) and Gj(0,T) indicate that each bilayer has a well defined BCS-like DOS, 
which coexists with an anomalous finite DOS at the Fermi energy that persists in the zero 
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FIG. 3: (a) A is normalized by its zero temperature value Ao plotted vs. T normalized by T c . The 
solid line is the BCS form. The inset magnifies the high T/T c range, (b) 770 obtained from the fits 
in Fig. 2c and 77 obtained from Gj(V) as a function of T c . The dashed line is the BCS value for 77. 

temperature limit. 

The existence of a zero temperature DOS at Ep in a homogeneous s-wave superconduc- 
tor without time reversal symmetry breaking perturbations is unexpected. In principle, it 
could arise for bilayers outside of the Cooper limit, for which is thick enough that the 



superconducting order parameter varies perpendicular to the bi 
depression of the bilay ers' T c 's with normal layer thickness j?l 



ayer plane. The exponential 
14j as well as previous work 



o. Ph/Cu systems 0, however, g ive us eo^ee th at the Ma yers studM he, are in t he 
Cooper limit. 

The "super weak coupling" behavior becomes evident upon closer inspection of the A(T). 
First, a scaled plot, A(T)/A versus T/T c reveals that the data for the three lowest T c 
bilayers fall below the BCS universal curve near T c . These deviations, which become larger 
for lower T c , are opposite those exhibited by strongly coupled superconductors. Second, the 
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gap ratio, 77, decreases from above 3.52 at the highest T c , to well below for the lowest T c 
bilayers. This behavior is shown in Fig. 3b where rj has been calculated using both Ao, t/o 
and A(0), 77 The difference between the two reflects an uncertainty inherent in defining 77 



for a system with a distribution of energy gaps. Each yields a different "averaged" value of n . 
Uncertainties in 77 due to T c determination are negligible for the very sharp transitions j?!. Il4j . 
The relatively good agreement between the measured A(T) and the BCS form, which yields 
the extrapolated A , however, leads us to believe that r/ provides the more reliable average. 
Whichever the case, it is clear that 77 falls well below the BCS weakly coupled value at the 
lowest T c 's. 

Super wea Hy C0Upled is unexpeetedQ a, th0Ugh Unto of it have b ee n 

observed previously. Merchant and coworkers reported a steady decrease in 77 to just below 
3.52 in their lowest T c SN bilayersQ] based on extrapolations of A from 1.5 K to zero 
temperature presuming the BCS temperature dependence. Similarly, 77 systematically falls 
with increasing sheet resistance in amorphous Bi films near the superconductor to insulator 
transition In addition, 77 <3.52 has been reported for alloys such as Af& 3 SVi[19], but the 
superconducting properties of these systems have been shown to be very disorder dependent 
and the reduced gap ratio values were most likely due to a damaged layer at the tunnel 
junction interface 201 ]. 

We suggest that the appearance of super weak coupling in SN bilayers is related to the 
subgap DOS. Within the current model of the subgap states, Gj(0,0) is proportional to 
the density of quasiparticles that are localized within the N regions [3, UJJ, lllj . These 
states cannot contribute to pairing at low temperatures. They can, however, contribute to 
pairing at high temperatures where the divergence of the coherence length allows them to 
become untrapped. Thus, the effective DOS available for pairing is higher near T c than at 
low temperatures. This effect can reduce 77 since A and T c depend on the low and high 
temperature DOS', respectively. 

We can test the above scenario semi-quantitatively by estimating the difference in the 
high and low temperature DOS implied by the measured 77. We use the Cooper limit relation 
for T c Q: 

where A is the electron-phonon coupling constant and (5 is a parameter that depends on the 



S 



0.20 



0.15 



0.10 



0.05 



0.00 



□ 

D ♦ 



♦ 5n 



□ ♦ 



□6/(0,0) 



□ 



□ 



1 2 
T c (K) 



FIG. 4: Calculated fractional reduction in the low temperature DOS, 5n, and the extrapolated 
zero bias conductance versus T c . 

Fermi velocities of the N and S layers and the interlayer coupling We take A to have 
a similar form, but with a coupling constant, Aa, that depends on d^. Since the A's are 
proportional to the DOS available for pairing, the fractional difference in the high and low 
temperature DOS' is given by Sn = 1 — Aa/A. Thus, r\ and Sn are related by: 

Sn = 1 - ( 1 - f sX log{r]hbu,ik)\ 
\ (3d N + d s J 

where n^f. is the gap ratio for which At = Aa- Sn is compared with G,(0,0) in Fig. 4. It 
was calculated using rj from Fig. 3c , the bulk Pb values rj^ik = 4.38 and A = 0.57, and 
/3/X — 0.4 obtained from a fit to T c (d^)[7\. The growth in the DOS at Ep matches well 
with the predicted reduction in the density of states available for pairing and thus, supports 
the scenario that quasiparticle trapping leads to super weakly coupled superconductivity in 
these bilayers. 

The data on ultrathin Pb-Ag bilayers presented here reveal a superconducting ground 
state with n decreasing below the BCS universal value as a finite DOS at the Fermi energy 
grows. Our analysis suggests that these features are related and provides further support 
for the picture that a growing fraction of quasiparticles decouple from the superconducting 
phase with decreasin g bi layer T c j^. Interestingly, 77's well below 3.5 and states within the 



superconducting gap 



18, 



2l| appear with decreasing T c in homogeneous thin films near the 



Superconductor to Insulator transition. The spontaneous formation of non-superconducting 



regions embedded in superconducting regions |22j, |23|, |24j may lead to quasiparticle local- 
ization and states at the Fermi energy in these pure, but highly disordered systems. Taken 
altogether, experiments suggest that super weakly coupled behavior and quasiparticle de- 
coupling may be general features associated with the destruction of these superconducting 
states. 
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